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ABSTRACT 

Context. The isocyanic acid (HNCO) presents an extended distribution in the centers of the Milky Way and the spiral 
galaxy IC342. Based on the morphology of the emission and the HNCO abundance with respect to H2, several authors 
made the hypothesis that HNCO could be a good tracer of interstellar shocks. 

Aims. Here we test this hypothesis by observing a well-known Galactic source where the chemistry is dominated by 
shocks. 

Methods. We have observed several transitions of HNCO towards L1157-mm and two positions (Bl and B2) in the blue 
lobe of the molecular outflow. 

Results. The HNCO line profiles exhibit the same characteristics of other well-known shock tracers like CH3OH, H2CO, 
SO or SO2. HNCO, together with SO2 and OCS, are the only three molecules detected so far whose emission is much 
more intense in B2 than in Bl, making these species valuable probes of chemical differences along the outflow. The 
HNCO abundance with respect to H2 is 0.4-1.8 10~* in Bl and 0.3-1 10~^ in B2. These abundances are the highest ever 
measured, and imply an increment with respect to L1157-mm of a factor up to 83, demonstrating that this molecule is 
actually a good shock tracer. 

Conclusions. Our results probe that shocks can actually produce the HNCO abundance measured in galactic nuclei and 
even higher ones. We propose that the gas phase abundance of HNCO is due both to grain mantles erosion by the 
shock waves and by neutral- neutral reactions in gas phase involving CN and O2. The observed anticorrelation of CN 
and HNCO fluxes supports this scenario. The observed similarities of the HNCO emission and the sulfured molecules 
may arise due to formation pathways involving also O2. 

Key words. ISM: individual (L1157) - ISM: jets and outflows - ISM: molecules - stars: formation - shock waves 
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1. Introduction 

Interstellar isocyanic acid (H NCO) was first detected to- 
wards Sg r B2 (Snvder & Buhl l972HChurchwell et al.lll98()l: 
iKuan fc Snvderi I1996f) . Since the first detection in this 
source, the molecule has been observed in other hot 
cores around massive (jBlake et al.lll987t iMacDonald et al 
1996 ) and low ma ss protostars (|van Dishoeck et al 



mapped: IC342 



I995I : iBisschop et al.l I2008D. It has also b een detected 
in translucent clouds fTurner et al.l 119991) and in the 
dense regions of Galact i c mol ecular clouds packson et al.l 



119841: IZinchenko et all l2000l), inc l uding those in th e 
Galacti c cente r (iHiittemeister 1993; Lindavist ct al. 199J 

Rizzo ct al. 2000; Minh ct al. 2001 

20081 ). HNCO has also been detected 



jPah mcn et a.1 

' iMartm et all ^ 

some extragalactic sources (iNguven-Q-Rieu et al.l 119911 ; 



iMeier fc Tumeil 120051; IMartm et all l2009f) . The isotopo- 
loguc H CNO (fu l minic acid) has recently been detected by 
iMarcclino et all (|2009D . 

The emission in HNCO lines with iC.i = OQ is clearly 
extended in the two Galactic nuclei where it has been 
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^ HNCO is a planar, nearly linear, slightly asymmetric prolate 
rotor. The notation for a given level is Jk_iKi ■ The K-i = are 
usually excited thermally, while the excited K ladders (K-i > 
0) ar e probably excited by FIR radiation l|Churchwell et ahl 
119861 ). 



(jLindavist et al. 



Meier fc T urner"2 005D and the Milky Way 
1995HDah mcn ct aill997t) . In the Galactic 



center the spatial distribution of the HNCO emission is dif- 
ferent from that of most other tracers. In particular, there is 
a HNCO peak witho ut an associated CO peak at Galactic 
longitude / = 1.65° (jPahmen et al.l[l997l) . This very spe- 
cial distribution suggests that HNCO can be an important 
tracer of some physical processes that are not well revealed 
by other molecules. It is possible that the molecule is trac- 
ing shocks, which are thought to take place at this Galactic 
longitude due to the gas dynamics in the barred poten- 
tial o f the Milky Way (Rodriguez- F ernandez fc CombesI 
I2OO8D . Indeed, in the Galactic center, the highest gas 
phase abundances of SiO (a well-known shock tracer 
are measured in this region ( Hiittemeister et al.l |199< 



I egic ^ 

iRodriguez-Fernandez et al.ll2006[ ). Also on the basis of a 
special morph o logy of the HNCO emission in Sgr B2, 
iMinh fc Irving (|2006f ) suggested that HNCO is enhanced 
by shocks. Regarding IC342, the HNCO emission resem- 
bles that of CH3OH. In particular, HNCO is detected 
not only in the nuclear ring but also in the inner spi- 
ral arms or dustlanes. Since the dustlanes are thought 
to be the l ocus of strong shocks in barred gala. xies (see 
for instance [Rodriguez- Fernandez fc Combesll2008l and ref- 
erences therein). iMeier fc Turneij ( 2005() suggested that 
HNCO, as CH3OH, could trace large scale shocks. 
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Fig. 1. Integrated C0(2- l) emission from the L1157 out- 
flow (jBachiller et al.ll2001[ ). We show the position of L1157- 
mm (star), Bl (empty square) and B2 (filled square). Bl 
and B2 are the positions of the sh ocks as traced, for in - 
stance, by SiO(3-2) (see Fig. 6 of iBachiller et alJl200lD . 
The CO emission peaks behind the actual positions of the 
shocks. The different circles represent the half maximum 
contours of the IRAM 30m beam at the frequencies of the 
lines discussed in this paper (see Tabled]). 
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Table 1. Observational parameters. 



Line 




Be//'' 




FWHM" 


-L sys 


rms 




GHz 






arcsec 


K 


mK 


4-3-' 


87.925 


0.78 


0.95 


28.4 


100 


12-18 


5-4" 


109.905 


0.75 


0.95 


22.3 


180 


16-19 


6-5'^ 


131.885 


0.69 


0.93 


18.8 


215 


28-34 



Notes. Frequency Beam efficiency Forward efficiency 
''^^ Full width at half maximum of the telescope main beam 
System temperature 4o,4-3o,3 5o,5-4o,4 6o,6-5o,5 



2. Observations 

We have observed three lines of the K^i — ladder of 
HNCO towards the two main shocks in the southern lobe 
of the L1157 molecular outflow: the Bl and B2 positions 
of'Bachi ller &: Perez G utierrez (1997). In addition, we have 
observed as reference the continuum source L1157-mm (the 
protostar). The equatorial coordinates of L1157-mm are 
RA=20^39'"06.19" Dec=68°02'15.9", (J2000). The offsets 
of Bl and B2 with respect to L1157-mm are (20", -60") 
and (35",— 95"), respectively. 

The observations were done with the IRAM 30m tele- 
scope in Pico Veleta in July 2007. The line quantum num- 
bers, frequencies and the telescope parameters are listed in 
Table [TJ As backends we used the 100 KHz filterbank and 
the VESPA autocorrelator with a channel resolution of 20- 
40 MHz, which allows to study in detail the line wings. 
The observations were carried out in position switching 
mode with the OFF position located 24' NE from L1157- 
mm. Typical system temperatures and rms noise of the 
spectra are are given in Table [TJ 

Table [5] gives the integrated flux of the different lines 
computed from spectra in T,nb units (forward and beam ef- 
flciencics arc given in Table[T]). The uncertainties in the cal- 
ibration of the IRAM 30m telescope ar e within 10 % at the 
frequ ency of the lines discussed here (jMauersberger et al.l 
Il989h . 



3. Results 



Nevertheless, the hypothesis that HNCO is a good shock 
tracer at the scale of galaxies still needs to be probed since 
the HNCO emission has never been studied in well-known 
Galactic templates of interstellar shocks. In order to bet- 
ter understand the excitation and the chemistry of this 
promising molecule, we have observed the protostar LI 157 
and its associated molecular outflow. T his outflow presents 
the morphological s ignature of shocks (Gueth et al.lll99"8t 
ICodella et al.l l2009t ) and it is fre quently used to ben ch- 
mark numerical models of shocks (jCusdorf et al.l[2008al lb[). 
L1157 is the best example of "chemically active outflow" 
(jBachiller et al. 2001) and a template of shock chemistry, 
since many species exhibit l arge abundance increments 
with respect to the protostar (jBachiller fc Perez Gutierrei 
[T997I) . Therefore, it is a source of choice to characterize 
the emission of a given molecule in a shocked environ- 
ment (see for instance Bachiller fc Perez Gutierre3 1997^ 



Bachillcr ct al. 2001; Bc nedettini et al 
2008c ,Codella et al..,2009i ). 



20071: lArce et alj 



3.1. Line profiles: differences between the observed positions 

Figure H) shows the spectra. The 4o4 — 3o3 and 5o5 — 4o4 
lines have been detected in all the positions. In addition the 
606 — 5o5 line has also been detected toward L 1157mm. The 
lines are narrow towards L1157-mm (linewidths of 0.7-0.9 
km s~^) but the line width increases by more than a factor 
of 3 towards the Bl and B2 positions. In particular, they 
show a prominent blue- wing towards Bl. 

Figure [3] shows a comparison of the HNCO (5o5 — 
4o4) line profile measured towards the three observed po- 
sitions and the spectra of other molecules observed by 
IBachiller fc Perez Gutierrez! (jl997f ) . The line profiles are in 
overall good agreement with the other molecules. In par- 
ticular, the HNCO profiles are very similar to those of 
CH3OH, H2CO, SO and SO2, which are all shock tracers 
whose abundance has increased by more than a factor of 
10 in the shocked gas with respect to the quiescent gas in 
L1157-mm. The HNCO line wings are na rrower than SiO 
lines, in agreement with the findings of IZinchenko et al.l 
(j2000t ) in Galactic dense cores, but the SiO line profiles 
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are specially broad in the L1157 molecular outflow, being 
broader than those of all the other molecules as well. 

As discussed by iBachiller fc Perez GutierreS (|1997D . 
most of the molecules exhibit lines with similar intensity 
in B2 and in Bl. The exceptions are CN, which is very in- 
tense towards Bl but it is not detected towards B2, and 
the sulfur-bearing molecules SO2 and OCS, whose lines are 
much more intense in B2 than in Bl. It is interesting to 
point out that HNCO, whose lines are a factor of 2 more 
intense in B2 than in Bl, is indeed the only non sulfur- 
bearing molecule showing this behavior. 

3.2. Enhanced HNCO abundances in the shocked gas 

The critical density of the HNCO transitions discussed here 
is high ( 5 10^-10^ cm~^; derived using the collisional coef- 
ficients in the LAMDA databas e for temperatures from 20 K 
to 320 K. ISchoier et al.ll2005[) . Therefore, it is likely that 
most of the emission comes from dense and small clumps 
as those foun d in Bl with the TRAM Plateau de Bure in - 
terferometer (jBenedettini et al.ll2007l : ICodella et al.ll2009t) . 
For instance, the angular size of the CH3OH clumps is 10". 
However, without high angular resolution observations it 
is difficult to assign a very precise value to the size of the 
source emitting in HNCO. Therefore, all the following re- 
sults have been obtained in two limiting cases: first, assum- 
ing that the size of the emitting region is 10" and computing 
source brightness temperatures taking into account the cor- 
responding beam filling factor. Second, assuming that the 
emitting region fills the primary beam of the telescope, and 
thus using Tmi, temperatures. 

3.2.1. LTE analysis 

We have computed HNCO column densities assuming opti- 
cally thin emission and LTE excitation. We have computed 
excitation temperatures from the 5o5 and 4o4 populations 
(Tg^) and from the 606 and 4o4 populations (T^;^). The tem- 
peratures Tg^ are very low (5-6 K for mm and < 6 K for 
Bl and B2). The excitation temperature T^^^ varies from 8 
to 11 K for L1157-mm (depending on the assumptions on 
the source size), while it varies from 9 to 19 K in Bl and 
B2 (see Table [31). In order to study with great precision 
the HNCO excitation more than three transitions would be 
needed. Here, in order to obtain a conservative estimation 
of the total column density of HNCO we have extrapolated 
the population in the 5o5 and 4o4 levels using the T^;^. Of 
course, the total HNCO column density would be higher 
if the excitation temperature is actually lower. Using T^j^ 
and assuming LTE excitation, the total HNCO column den- 
sities are 0.2-0.8, 1-1.9 and 1.5-3.8 (in units of 10^^ cm~2) 
for mm, Bl, and B2, respectively. Since the excitation is 
similar in Bl and B2, the highest line intensities measured 
towards B2 translate into larger column densities, in con- 
trast t o the CO column density which is a factor of 2 lower 
in B2 (jBachiller fc Perez Gutierredll997D . 

3.2.2. RADEX analysis 

In addi tion, we have studied th e HNCO excitation using 
RADEX (jvan der Tak et al.ll2007t) , which is a non- LTE exci- 
tation and radiative transfer code that decouples the sta- 
tistical equilibrium and the radiative transfer equations us- 



ing the escape probability method. The calculations assume 
lines of a given width and rectangular shape. Therefore, the 
change of the optical depth over the profile is not taken into 
account. This is not a problem if the opacity of the lines is 
not very high. We have used the HNC O-H2 collisional coef - 
ficients listed in the LAMDA database ()Sch5ier et al. I [20051 ). 
The computations use a background temperature of 2.73 K. 
We have modeled the line fluxes (or integrated intensities) 
as a function of the hydrogen volume density (71^2)1 the 
HNCO column density (Nhnco), the kinetic temperature 
{Tk) and the velocity dispersion. The range of hydrogen 
densities used is lO'^ — 10^ cm~^ and that of Nhnco is 
10^^ — lO'^^ cm~^. Two kinetic temperatures have been 
considered: a low kinetic temperature (30 K), only slightly 
higher than the excitation temperature derived with the 
HNCO lines and a much higher temperature (250 K), since 
the actual gas temperature in the shocked gas can be con- 
siderably higher than the HNCO excitation temperatures. 
These temperatures probably cover the different physical 
conditions of L1157-mm and the shocked sources Bl and 
B2. We have computed models with velocity dispersions 
(linewidths) from 1 to 6 km s~^, which cover the mea- 
sured range of linewidths in the three observed positions. 
In the range of physical conditions compatible with the 
observations, the opacity at the peak of the lines is low. 
Therefore, the predicted integrated intensities are indepen- 
dent of the actual linewidth used to compute the models. 
For instance, the figures discussed below have been com- 
puted with a linewidth of 3 km s^^. 

We have based our analysis on the two lowest lines, 
which are detected towards the three observed positions. 
Figure SI shows the model predictions for the 5o5 — 4o4 to 
4o4 — 3o3 line ratio and the 4o4 — 3o3 brightness temperature 
(for a source size of 10") as a function of uh^ and Nhnco- 
The same diagrams for a source size equal to the telescope 
main beam (using Tmb) are shown in Fig. [SI 

The kinetic temperature cannot be constrained since 
equally good solutions can be found just changing the hy- 
drogen density by 0.5 dex. For a given source size and 
kinetic temperature, the H2 density is similar for Bl and 
B2. The density varies from lO'' — 10^'^ cm~'^ for a source 
size of 10" and Tk = 250 K, to IQ'^-^ - lO^ ^ cm^^ for 
sources filling the main beam and Tk — 30 K. In L1157- 
mm, is lower than that in Bl and B2 by 0.2 — 0.3 dex. 

The total HNCO columns densities are 0.2-1.3, 0.6- 
2.5 and 1.3-5.0 (in units of 10^"^ cm~^) for mm, Bl, 
and B2, respectively (Table [Sj. The lower limits are sim- 
ilar to those derived using LTE, while the upper lim- 
its for Bl and B2 are somewhat higher because opacity 
effects start playing a role. Table [3l also shows HNCO 
abundances computed with the RADEX estimations of the 
HNCO column de nsity and the total H2 column densi ty de- 
rived from CO bv IBachiller fc Perez Gutierrez! (|1997f ). The 
HNCO abundances are (0.3 - 1.2) x lO'^ for L1157-mm, 
(4.3 - 17.9) X 10-^ for L1157-B1 and (25 - 96) x 10"^ for 
L1157-B2 (Table 121). Therefore, there is a clear increase of 
the HNCO abundance in the shocked regions, by a factor 
of 6-14 in Bl and by a factor of 34-83 in B2. This is the 
first direct evidence that HNCO is indeed a good tracer of 
shocked gas. 

In the range of physical conditions and HNCO column 
densities derived from the 5o5 — 4o4 and 4o4 — 3o3 lines, the 
model predictions for the 606 — 5o5 to 5o5 — 4o4 ratio is 0.6- 
0.7. In contrast, the measured ratio is 0.3-0.4 for L1157-mm, 
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Fig. 2. Observed HNCO spectra. From bottom to top the different panels show the HNCO 4o4 — 3o3, 5o5 — 4o4 and 
6o6 — 5o5 spectra towards L1157 mm, Bl and B2. 



and it is lower than 0.3 for Bl and B2. As already mentioned 
above when discussing the LTE analysis, a multi-transition 
survey would be needed to fully understand the excitation 
of the molecule. In addition, HNCO maps with high angular 
resolution will also be needed to determine the positions 
and the sizes of the different HNCO clumps, in particular 
to see if there could be some HNCO emission out of the 
IRAM-30m beam at 2mm. 



3.2.3. Comparison with other Galactic sources 

Table 2] shows HNCO abundances with respect to H2 mea- 
sured other Galactic sources as translucent clouds, photon- 
dominated regions (PDRs), dense cores and hot cores. The 
highest HNCO abundance is found in hot cores and dense 
cores: up to - 8.7 IQ-^. The HNCO abundance in L1157- 
mm is similar to that measured in those dense cores and 
hot cores with the lowest HNCO abundance and in the only 
"hot corino" where the molecule has been p reviously de- 
tected, IRAS 16293 (Ivan Dishoeck et al.lll99,'ih . The HNCO 
abundance in the L1157-B1 shock is comparable, but can 
even be higher than that in hot cores and dense cores. 
Regarding the L1157-B2 shock, the HNCO abundance is 
at least 4-11 times higher than that in hot cores and dense 
cores. The HNCO abundance in the L1157-B2 shock are 
indeed the highest HNCO abundances ever measured. 

4. Discussion 

4.1. HNCO formation in molecular clouds and hot cores 

The possible formation pathways of HNCO in shocks have 
never been modeled. In contrast, HNCO has been included 



Table 2. Measured fluxes 



Line 




mm" 


Bl" 


B2'= 






K kms-^ 


K kms^^ 


K km 


4o.4 - 


3o,3 


0.19 


0.71 


1.24 


5o,5 - 


4o,4 


0.19 


0.84 


1.38 


60,6 - 


5o,5 


0.073 


< 0.25 


< 0.20 



Notes. Computed in the velocity range from: 1.1 to 3.7 
kms"\ -4.8 to 5.0 kms"\ -2.6 to 6.1 kms'\ 



Table 3. Resuhs 



Source 


r-p54 a 
J- ex 

K 


Nhjvco'' 
10^^ cm"^ 


10" cm"^ 


X(HNCO)'' 
10-s 


mm 


8-11 


0.2-1.3 


0.2-1.3 


0.3-1.2 


Bl 


11-19 


0.6-2.5 


0.6-2.5 


4.3-17.9 


B2 


9-16 


1.3-5.0 


1.3-5.0 


25-96 



Notes. "-"^ Excitation temperatures derived from the popula- 
tions in the 5o,5 and 4o,4 levels. HNCO column density as- 
suming LTE, optically thin lines and extrapolating using T^^. 
'•"^ HNCO column density computed using radex. HNCO 
abundan ce with respect to H2 (using the H2 column densities 
given bv lBachiller fc Perez Gutierrezlll997l ) 



in some model s of dark cloud s and hot core chemistry. In 
the models bv llglesiai (|1977[ ). HNCO is produced by the 
ion-neutral reaction H2 + HNCO+ H2NCO+ + H fol- 
lowedby H2NC0+-He- ^ HNCO-^H. lTurner et al](ll999ll 
have considered three possible formation pathways, among 
them the only efficient one is the neutral-neutral reaction 
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CN + O2 ^ NCO + followed by NCO + H2 ^ HNCO + H. 
The last reaction lias an activation barrier that has been 
estimated to be ^ 1160 K byjurner et al. (1999,) but that 
could be as high as 4465 K (fTideswell et alJ 12010 ) . In any 
case, this reaction is not efficient at the typical tempera- 
ture of a hot core (~ 200 K). Therefore, gas phase chem- 
istry alone cannot explain the HNCO abundances in hot 
cores. This is also the case for other species as the com- 
plex organic molecules. Chemical models have consequently 
been developed to include reaction s on the dust grains sur- 
face jCaselliet al . 1993; Garrod et a l. 2008: Tidesweh et al.l 
|2010() . In those works, hot cores are modeled in two stages. 

In the first phase, a dark cloud of uh ^ 10"* cm~^ 
suffers an isothermal collapse. The collapse phase is halted 
once uh ~ lO'' cm~'^, which occurs after approximately 
5 10^ years (free-fall time). The typical temperature at this 
phase is 10-20 K and surface chemi stry is very impor- 
tant then. Surface ch emistry models ()Garrod et al.l l2008t 
iTidesweh et"alll2010D produce HNCO as a secondary rad- 
ical formed via the reaction NH CO HNCO.Q On 
the grain surface, HNCO can r each a maximum abundance 
of 10^^ (jTideswell et al.l I2OIOD before being destroyed by 
new reactions with primary radicals (H, CH3, HCO, NH, 
...) giving c omplex species as HNCHO, HNCOCHO, or 
CH3CONH (iGarrod et al.ll2008D . 

The cloud-collapse phase is followed by a warm-up 
phase in which the temperature increases to ^ 200 K. 
The timescale of the warm-up phase can be as short as 
~ 5 10^ yr in hot cores but it can reach 10^ yr in their 
low mass equivalent (hot cor inos), giving a simil ar although 
slightly different chemistry (jCarrod et al.ll2008h . 

During the warm-up phase molecules are evaporated 
from the grain mantles. W ith respect to HNCO, the mod- 
els bv lGarrod eT^ (|2008D and lTidesweh et al.l (|2010D show 
that its abundance increases monotonously with time in the 
gas phase, which implies that HNCO is not directly ejected 
from the grains. Instead, HNCO is formed by the destruc- 
tion in the gas phase of the comple x molecules formed 
from HNCO on the grain surface (|Garrod et al.l l2008t 
I Tides well et al.l |2010D . This mechanism can explain the 
HNCO abundances measured in hot cores of 10~^ — 10~^. 
As alr eady mentioned, th e contribution of the lTurner et al.l 
()1999() and llglesiasi (|1977D formation pathways is negligible 
at the moderate temperature of a hot core. 

The abundance of HNCO in L1157-mm is similar to 
that measured in the hot corino of the only other low mass 
protostar where HNCO has been observed (IRAS16293- 
2422, see Table g]). Therefore, the HNCO abundance mea- 
sured in L1157-mm can also be explained in the context 
of the hot-cores models discussed above. This would imply 
that L1157-mm can be considered as a hot corino, which 
is indeed in agreement with the intense emission of water 
detected by Herschel towards mm (Nisini et al. in prep). 
However, more observations will be needed to clearly es- 
tablish the hot corino character of L1157-mm, in particular 
observations of complex organic molecules. 



^ HNCO formation on the grain surfaces could also occur 
by reactions of OCN^ with NHJ or HsO^. The inverse reac- 
tions (HNCO with NH3 or H2O) are inv oked to explain the 
presence of OCN~ ices in th e dust grains, l|Demvk et al.lll998l : 
Ivan Broekhuizen et al.ll200l ). To our knowledge, these reactions 
have never been included in chemical models. 



4.2. HNCO formation in shocks 

What is the formation pathway of HNCO in the L1157 
outflow shocks? The most likely scenario is a combina- 
tion of grain surface and gas phase chemistry. In the con- 
text of shocks, molecules formed on the grain surfaces can 
be ejected to the gas phase due to grain sputtering in- 
stead of thermal evaporation as in hot cores. Therefore, 
the chemistry will be sensitive to the grain mantles com- 
position at the time of the shock arrival. Currently, there 
are no models that study the HNCO abundance in dark 
clouds including grain s urface chemistry. Dark clouds mod- 
els by iTideswell et al.l (|2010f ) only consider gas phase re- 
actions. Therefore, the exact grain mantle composition be- 
fore the arrival of the shock in L1157-B1 and L1157-B2 is 
not known. Nevertheless, one can compare with the cloud- 
collapse phase of hot core models (the main difference is 
that the density during the collapse phase reaches higher 
values than in dark clouds). As discussed in the previous 
section, during the collapse phase the HNCO abundance 
on the grain surfaces has a peak of 10~^, grain sputtering 
at that time will immediately give rise to very high HNCO 
abundance in gas phase. Even the very high HNCO abun- 
dance in L1157-B2 (^ 10^^) could be accounted for in this 
context. 

On the other hand, gas phase chemistry in shocks dif- 
fers considerably from that in hot cores because the tem- 
perature in the shocked gas can be much higher than that 
i n a h ot core. As already pointed out by IZinchenko et al.l 
(j2000[ ). in shocks the reaction NCO -I- H2 can be efficient 
in spite of the activation barrier. In addition, the reaction 
CN O2 ^ NCO + O wiU be favored by the h igh O2 abun- 
danc es that are pred icted in post-shock gas (jBergin et al.l 
IT99I iGusdorf et al.ll2008allbl) .Pi 

Using ITideswell et al.l (|2010f ) results, it is possible to 
verify if gas phase chemistry alone could explain the HNCO 
abundance measured in the L1157 shocks. In one of their 
dark cloud models (DCl), they have studied the HNCO 
formatio n exclusively via the gas phase neutral-neutral re- 
actions of lTurner et al.l (jl999t ) as a function of the rate coef- 
ficient for the reaction NCO-KH2 -J> HNCO + H. Their Fig. 
1 shows that a rate coefficient k higher than 10"^** cm'^ s~^ 
is needed to explain the observed abundances in LI 157 of 
10^^ — 10^^. However, to have high HNCO abundance (a 
few 10^*) in a short time (several 10'^ yr) the rate coefficient 
must be higher than 10~^^ cm^ s~^. Taking into account the 
formula k(T) = 1.5 10 ~" exp(-4465/T) cm'^ s'^ given by 
ITideswell et all (|2010D . k = 10"^^ cm^s"^ implies a tem- 
perature higher than 1600 K. In contrast, the minimum 
temperature would be 428 K if t he activation barrier is 
only 1160 K (iTurner et al.l[T99l . IGusdorf et aP (l2008bD 
have recently computed shock models to compare the mod- 
elled SiO and H2 emission to observations of LI 157. The 
best agreement between models and observations is found 
for a preshock density of 10** cm~^ and a shock speed of 
20 kms~^. In such a shock, the gas temperature is higher 
than 430-1600 K for a few thousand years, which is con- 



^ This does not mean that the current O2 abundance in the 
L 1157 shocks is high, but it may have been high just after the 
passage of a shock-wave as predicted by shock models. Since the 
chemistry in shocked regions is fast, O2 can be rapidly converted 
into other species. Therefore, there is not any contradiction with 
the fact that the O2 abundan ce in the interstellar medium has 
been found to be low (< 10"'^ iLarsson et al.ll2007r i. 



5 



N. J. Rodriguez-Fernandez et al.: HNCO enhancement by shocks in the L1157 molecular outflow 



Table 4. HNCO abundance with respect to H2 in different 
astronomical environments. 



Source 



X(HNCO) 



References 



L1157-B2 
Translucent clouds 
Orion bar PDR 
Dense cores 
Orion Hot core 
SgrB2M 
W3(H20) 
G34.3+0.15 
IRAS 16293 
GC clouds 
Starburst galaxies 



25-96 


1 


0.2-5 


2 


< 10-^ 


3 


0.2-8.7 


4 


5 


5 


1 


6 


5.0-27 


7 


1.4 


8 


0.17-4.0 


9, 10 


3-48 


11,12 


1-6.3 


13 



Refer ences. (1) T his wo rk; (2) I Turner et al 
1999f) :f3) iJansen et al"] (I1995D : ( 4) iZinchenko et al 
20001): (5) iBlake et al.1 (1198711 : (6) iMinhet ajJ (Il998l) 



7) iHeimich fc van Dishoeckl l(l99W ): (8) [Bockclcc- Morvan et al.'l 
(2000), using th e HN CO column density given by 
iMacDo nald et al.l (119961) and the H2 column density 
given bv iHatcheU et al.1 (Il998f); (9 ) Ivan Dishocck et al.' 
19951): (10) iBisschop et all (I2008D : (11) Rizzo c t al.| 
200 0): fl2) iMartm et al.1 (|2008l ). and references within; 



l^l3) lMartm et al.l (|2009l )" 



sistent with the dynamical a ge of the Bl and B2 shocks 
(2000-3000 yr, iGueth et ai1 Fl998) and with the chemical 
time necessary to have an HNCO abundan ce of a few 10~^ 
only with gas phase reactions (Figure 1 of iTideswell etHI 
I2OICI ). Therefore, gas phase chemistry alone could explain 
the HNCO abundance in L1157-B1. However, it is unclear 
whether a HNCO abundance as high as that measured in 
L1157-B2 (~ 10~^) can be attained in a few thousand years 
only with gas phase reactions. 

Therefore, we reckon that the most likely explanation to 
the high HNCO abundances in L1157 is dust grain mantles 
processing by the shock waves followed by neutral-neutral 
reactions in gas phase. 

4.3. Chemical differences between Bl and B2 and the role of 
the O2 chemistry 

We have discussed the line profiles and the line intensities of 
different species in Sect. 13.11 The HNCO profiles are very 
similar to those of SO and SO2. In addition, the lines of 
most of the molecules show similar intensities in Bl and 
B2 with the exception of CN, which is more intense in Bl, 
and HNCO, SO2 and OCS, whose lines are more intense 
emission in B2 than in Bl. These similarities of HNCO and 
the sulfured molecules is somewhat puzzling. 

The actual reason of the chemical differences in Bl and 
B2 is not clear. The present gas density in Bl and B2 is 
simila.r (see Figs. P and [5] and iBachiller fc Perez Gutierred 
Il997t iNisini et al.ll2007[ ). In contrast, based on the line pro- 
files, the shock velocity can be higher in Bl than in B2. 
Taking into account the extreme SiO line wi ngs, the differ- 
ence co uld reach 10 km s~^, which after the lGusdorf et al.l 
(|20083) models c ould be significant. 

Alternatively, IBachille r fc Perez G utierred (jl997t ) have 
suggested that the chemical differences between Bl and 
B2 could be due to different shock ages. Indeed, in con- 
trast to SO2 and OCS, the H2S abundance is lower in 



B2 than in Bl. This is consistent since H2S is a parent 
molecule for other sulfured species. H2S is formed in the 
grain surfaces and released to the gas phase by effect of the 
shock waves. Other sulfur-bearing molecules like SO and 
SO2 are produced in gas phase v ery quickly (few 10'^ yr) 
via reactions with H, OH a nd O2 (jPineau des Forets et al.l 
I1993I : IWakelam et al.l 120051 ). One possible explanation of 
the differences in Bl and B2 is that the B2 shock could 
be older than that in Bl and that H2S has been converted 
into SO and SO2. Even if there are un certainties regarding 
the sulfur chemistry (see for instance ICodella et al.l l2005l 
and references therein), sulfur-bearing molecules like H2S, 
SO, SO2, H2CS and OCS have been invoke d to be poten- 
tial valuable probes of chemical evolution (ICodella et al.l 
I2OO5 I: 'Wakelam et al."2005t Berpin et al.l l2009[). In partic- 
ular, Wakelam ct al. (2003) have found that the SO2/SO 
ratio increases with the shock age. The higher SO2/SO 
ratio i n B2 compared to Bl also po i nts to an older shock 
in B2. IBachiller fc Perez Gutierred (jl997() have also pro- 
posed that the different CN abundance in Bl and B2 can 
be accounted for in this scenario of an older shock in B2, 
since after an enhancement of CN in the shock, reactions 
like CN-l-0 CO-I-N could be very efficient to destroy 
CN. Our data globally agree with this scenario, nevertheless 
the clear anticorrelation of the CN and HNCO intensities 
and the similar abundances of CN and HNCO in L1157- 
Bl suggest that all the CN can indeed be transformed into 
HNCO once the shock has increased the temperature and 
the neutral-neutral reactions form NCO from CN and O2 
and HNCO from NCO (see Sect.H^]). Therefore, the chem- 
ical differences from Bl and B2 do support the scenario 
proposed in Sect. 14.2] to explain the HNCO abundances in 
the L1157 shocks. 

In addition, the scheme proposed in Sect. 14.2! 
also gives insight on the possible link of HNCO and 
the sulfured molecules, in particular with SO and 
SO2. In shocks these molecules can be produced by 
the reactio ns S + O2 SO + O and SO OH -» • 

S Oo + OH (iPineau des Forets et al.l[T99l ICharnlevI [19971 : 



IWakelam et al.ll2005D . Therefore, the common link with 



HNCO would be formation pathways involving O2. 

Another important shock tracer that could be linked to 
the O2 chemistry is SiO. Nevertheless, the situation regard- 
ing this molecule is more complex. First, recent models can 
explain the SiO emission in shocks without a significant 
contribution of Si oxidation in gas phase, either by sput- 
tering in gas-g rain collisions if there is already SiO in the 
grain mantles ([Gusdorf et al.l 2008bl) or bv dust vaporiza- 
tion in grain-grain collisions ( Guillet et al.l [2009). Second, 
if one assumes that SiO is formed in gas phase, there is 
a threshold shock velocity of ~ 25 km s ~^ in order to 
eject S i from the grain cores by sputtering (jCusdorf et al.l 
i2008al ). In contrast, other molecules as the organics are in 
the grain mantles and there is not such a high velocity 
threshold for them to be ejected to gas pha se. Therefore, 
both if SiO comes directly from the grains (jCusdorf et al.l 
2008b: Guillet et al. 2009) or if it is formed in gas phase 
( Gusdorf et al.ll2008al ) . due to the shock velocity threshold. 



it is not expected to find similarities (and not observed) in 
the SiO emission and that of SO, SO2 and HNCO. 

It is interesting t o note that the im pressive SiO blue 
wing in Bl (Bachil ler fc Perez Gutierre z 199Z) suggests 
that SiO formation is indeed favored, in comparison to 
other molecules, in the higher velocity gas. In this con- 
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text, chemical differences in Bl and B2 could also be due 
to different shock velocities. If the shock velocity in B2 is 
actually lower than in Bl and the SiO formation less effi- 
cient, more O2 will remain available in B 2 to form HNCO, 
SO and SO2. 

4.4. On the origin of the HNCO emission in galactic nuclei 

The HNCO abundance in the molecular clouds in the 
center of the Milky Way and in the nuclei of starburst 
galaxies is similar or a bit higher (at most by a factor of 
2) than in Galactic hot cores (Table S]). HNCO becomes 
another piece in the well known puzzle of the chemistry 
of the Galactic center molecular clouds. This puzzle 
can be summarized as follows: the abundance of SiO 
and complex organic molec ules in the Galacic center 
is as high as in hot core s (iMartfn-Pintado et al. 19971 



Rodriguez-Fernandez et al. 1 120061: iReauena- Torres et al.l 



20061) . In contrast, (i) the emission in the Galactic center 
is extended over the central 300 pc and it does not 
resemble a collection of discrete sources with the size of 
a hot core (-- 0.1 pc), (m) the gas density (10^ - 10^ 
cm~^) is much lower than in hot cores (Hi) and the 
dust temperatures (< 30 K) are also much lower than 
in hot cores. The Galactic center clouds present a "hot 
core chemistry without hot cores" (|Reauena- Torres et al.l 
[2001 . The origin of this chemistry is not known although 
it is thought to be due to sort ie type of mechani- 
cal processes as shock waves (iMartfn-Pinta do et al. 
119971 I2OOII: iRodrfguez-Fernandez et al.l |2004 . The 
origin of the shocks can be related to the com- 
plex dynamics in t he in ne r regions of the Galaxy 
(Hiittcmcister et al.' 1998: Rodriguez- Fernandez et aD 
[2OQ6; Rodriguez-Fernandez & Combes 200^^ 

Based on the spatial distribution of th e HNCO and the 
compa rison with other species, as CH3OH, iMeier fc Tur^ 
(|2005[ ) have also suggested that the HNCO emission m 
IC342 is tracing shocks. Recently, iMartfn et all (|2008l 
[2OO&) have also proposed that shocks could be the expla- 
nation of the high HNCO abundances measured in galactic 
nuclei. 

Discussing the precise origin of shocks in Galactic nuclei 
is out of the scope of this paper. Nevertheless, our obser- 
vations support the scenario of HNCO tracing shocks in 
galactic nuclei since our LI 157 results probe a medium of 
moderate H2 density where the HNCO abundace is indeed 
high due to the grain processing and gas heating by shock 
waves. 



5. Summary 

We have observed three lines of HNCO towards the proto- 
star L1157 and its associated molecular outflow. Our aim is 
to characterize the HNCO properties in shocked gas. HNCO 
is well detected in the shocked gas, where the abundance 
increases by a factor of 6-34 with respect to the abundance 
in the protostar L1157-mm. The abundance in Bl and B2 is 
0.4-1.8 10~*^ and 0.3-1 10""^, respectively. The abundance 
in B2 is the highest ever measured, considerably higher 
than those in hot cores (a few 10~^) and galactic nuclei 
(10^^ - 10"*). Our results probe that the HNCO abun- 
dances measured in galactic nuclei can easily be attained in 
shocked gas, providing a solid basis to previous suggestions 



that the exte nded HNCO in galactic nuclei could trace lar ge 
scale shocks (Me ier k Turneill2005t iMinh k Irvinell2006[l . 

The dominant formation pathway of HNCO in hot cores 
is grain mantle evaporation of complex molecules formed 
from HNCO and subsequent dissociation to give again 
HNCO. In addition, there is contribution from gas phase 
reactions, but it is mi nor due to the h igh activation barri- 
ers of some reactions (jTideswell et al.| [2010). In the L1157 
shocks we propose a relatively similar formation pathway 
although with different relative importance of the dust sur- 
face and gas phase chemistry. First there should be grain 
mantles erosion by the shock waves that will increase the 
HNCO abundance in gas phase. Second, in contrast to hot 
cores, the neutral-neutral reactions in gas phase can be the 
dominant formation pathway due to the higher gas tem- 
peratures in the shocked gas with respect to the typical 
temperature of hot cores. These reactions start with the 
formation of OCN from CN and O2, which is expected to 
be very abundant in gas phase in a magneto-hydrodynamic 
shock of relatively low velocity (20-30 km s~^) like the 
L1157 shocks (Gusdorf et al. 2008a). The observed anticor- 
relation of CN and HNCO fluxes also support this scenario. 

HNCO line profiles are very similar to those of CH3OH 
and H2CO lines, which is not surprising since all those 
molecules have, at least partially, a common origin linked 
to the grain surface chemistry and subsequent desorption. 
In addition, HNCO line profiles are also similar to some 
sulfur-bearing molecules like SO and SO2. Another com- 
mon property of HNCO and the sulfur-bearing molecules 
is that they are the only species that exhibit a more in- 
tense emission in B2 than Bl. We propose that the common 
link of HNCO and the sulfured molecules is that in both 
cases there are formation pathways with reactions involving 
O2. HNCO together with the sulfured molecules are good 
probes of chemical differentiation in shocks. In particular, 
in the case of L1157 the highest abundances in B2 than in 
Bl could be due to a lower shock velocity and/or an older 
shock in B2 than in Bl. 
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Fig. 3. Comparison of the HNCO(5o5 — 4o4) line profiles (in 
red) with the spectra of other species published by Bachiller 
& Perez-Gutierrez (1997), in black hnes. The HNCO peak 
temperature has been scaled to the peak temperature of 
the other molecule in order to compare profiles. 
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Fig. 4. RADEX results for the 5o5 — 4o4 to 4o4 — 3o3 flux 
ratio (horizontal curves) and the 4o4 — 893 hne flux (verti- 
cal curves) computed with a brigthness temperature for a 
source size of 10" as a function of the hydrogen density, the 
HNCO column density and the kinetic temperature (30 K 
in the upper panel and 250 K in the lower panel) . A typical 
uncertainty of 10 % in the line flux has been taken into 
account. 
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